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ABSTRACT 


A technique is described for obtaining passive microwave measurements from space for remote 
sensing applications with hi^ spatial resolution. The technique involves measuring the product 
of the signal h'om pairs of antennas at many different antenna spacings, thereby mapping the 
correlation function of antenna voltage. The intensity of radiation at the source can be obtained 
from the Fourier transform of this correlation function. Theory will be presented to show how 
the technique can be applied to large extended sources such as the earth when observed from space. 
As an example, details vdll be presented for a system with uniformly spaced measurements. 
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I. INTRODUCTION 


Microwave radiometers in space possess a significant potential for monitoring the earth’s 
environment. Of current interest, for example, are microwave atmospheric sounders to complement 
infrared and visible sensors. Unlike infhued and visible radiation, microwave radiation can penetrate 
clouds and could provide meteorological data in the important regions beneath clouds and inside 
storms. However, microwave sensors in space require relatively large antenna systems and antenna 
^e ultimately imposes limits on the spatial resolution practical with such sensors. This paper 
d^cribes a technique to help overcome these limitations. 

The technique involves coherently measuring the product of the signal from pairs of antennas with 
many different antenna spacings. Each antenna pair measures the correlation function of the 
antenna volt^, R(x), at an argument, x, determined by the distance between the antennas. For 
distant sources, it can be shown that the correlation function is proportional to the Fourier trans* 
form of the intensity of the source at a frequency which depends on the spacing, x. Thus, by 
making measurements at many different spacinp one, in effect, determines the spectrum of the 
intensity of the source. The intensity can be obtained after aU measurements are complete by 
inverting the transform. The resolution obtained in this manner is determined by how well the 
COTielation function h^ been m^sured, not by the size of the antennas used. Thus, in principle, 
one can obtain very high resolution mai» of the source by measuring at m&.iy different baselines 
using relatively small antennas. 

This technique has been succemftiUy employed by radio astronomers to obtain very high resolution 
maps of radio sources in what is called earth rotation syntheris (Swenson and Mathur, 1 968; 

Brouw, 1973; Hewish, 1963). The Very Large Array m Socorro, New Mexico is an example. 
However, applying this technique to radimneters in space invohm some problems which are dif* 
ferait than those faced by the radio astronomer. For example, the objectt to be mapped in radio 
astfommiy are relatively small in anguMr extent and are quite distant whereas the earth is a‘ 


relatively large and nearby source for a radiometer in space, even one in geostationary orbit. 

It will be shown below that the problems of mapping nearby extended sources can be resolved by 
controlling the field of view of the radiometer by proper choice of the antennas and filters used 
in the correlator. Several versions of this technique have also been proposed for remote sensing 
from space. These include a proposal (Schanda, 1976, 1979) in which the baseline is changed by 
making measurements at different frequencies, a concept (Mel’nik, 1972) involving a moving 
system with a single fixed baseline in which m^urements must be made at several different time 
delays, and a modification which employs a scanning linear array which is rotated (C. Wiley, 
private communications). These systems can all be reduced to common principles which are the 
same as those employed in earth rotation synthesis. 

The purpose of this paper is to develop the mathematical foundation for the concept of aperture 
synthe^ as applied to remote sensing flt>m space. Several alternative configurations for its imple- 
mentation will be discussed, but the emphasis will be on the concepts and issues particular to 
microwave remote sensing from space rather than on the details of system design. In the sections to 
follow, tte response of a single antenna to radiation from the source win be computed first. Fhis 
win be done initiaUy in the frequency domain (Fourier transform of MaxweU’s equations with 
mpect to time) by employing the v^tor Helmholtz equation, tten the inverse transform wiU be 
taken to obtain the response in the time domain. In inverting the Fourier transform, it wiU be 
assumed that the bandwidth of the measurement is smaU compared to the nominal frequency of 
the system, an ^proximation which faciUtates the mathematics and is a realistic approximation 
for most microwave sensors. Next, this solution for the response of a sini^e antenna will be Used 
to obtain an expression for the spatial correlation function (le., the average of the product of the 
signato firom two antennas at dUfermt positions). In doing so, it will be assumed that the source 
and the maximum spacing between the antennas ate both small compared to the distance between 
the source and antennas, and it will be assumed that the source electric fields are incoherent. With 


these approximations, it will be shown that the intensity of the source can be obtained by taking 
a Fourier transform of the spatial correlation function. Methods for measuring the correlation 
function and meeting the restrictions imposed by the approximations made in the theory will be 
discussed. Finally, a hypothetical system in which the correlation function is sampled uniformly 
on a Cartesian grid will be analyzed and expressions will be d'/ived for the resolution and Held of 
view obtained with this sampling strategy. 
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II. ANTENNA RESPONSE IN THE FREQUENCY DOMAIN 

Imagine two parallel planes, the “scene” plane which contains a source of radiation and the “image’ 
plane where measurements are made. Also suppose an antenna is at position r in the image plane 
receiving radiation from the scene piane. In this section, it is desired to obtain an expression in 
the frequency domain for the output V(r, i>) from the antenna in terms of the electric field, 

Ejfr ', V ), in the scene plane. To begin, consider radiation frmn a small patch Ax' Ay' at r' in the 
scene plane. The electric field AE(r/r'; v) radiated from this small patch to the antenna can be 
obtained from the vector Helmholtz equation by assuming that the patch is a small aperture in 
an otherwise opaque screen on whkh the tangential components of the fields are zero. Following 
standard procedures (e.g., Tai, 1971), one obtains: 


AE(r/f';i;) 


[n X V X Ejlr', v )] • G(r/r ') + 
n xE^r',v)3‘VxG(f/f)| Ax' Ay' 


( 1 ) 


where » z is the unit vector normal to the scene plane and Gfr/r') is a dyadic Green’s function 
satisfying the free space wave equation in the region z > 0. A convenient choice for S(r/f ') is: 

G(r/r') ■ [go(r/r')-go(q/f')] I (2) 

where I is the unit dyadic, go(r/r') * e^^^^ ^4 it IF • r'l Is the free space scalar Green’s function 

and7i«r»2z^isthe “image” of the antenna’s position behind the scene plane. Now, the re* 
9onse of an anteima to radiation from the small patch can be expressed in terms of the voltage 
transfer (Unction of the antenna, A (r/f v ), as follows: 


A V(r/f ', p ) « AE(f/r', » ) • A(r/f ', » ) 


■ (3) 


The volt^ transfer function is a vector quantity to amount for the polarization of the antenna 
(C<^ and Zucker, 1969) and can be obtained from the reciprocity theorem for antmuias: One can 
show that AfF/f*, » ) ■ A*(r/r v ) ■ A^ Pg(r/r> where A« is the effective area of the receiving 
antenna and is the normalized power pattern of the antenna (e.g., Kraus, 1966). 
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To obrain the re^oiue of the antenna to all sources, it is now necessary to sum over all of the 
scene plane. In the limit of inflnitesimally small patches the sum approaches an integral and one 
obtains: 

V(f, V) ■ 

scene 


II 


AE(r/r', v ) • A(r/f', v) dr ’ 


ORIGINAL 
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■ JJ " x?xEs(f'i?)] G(f/f) + 
scene 

♦ [ n X Ej(r e )] • ? X Gfr/r*) | • A(r/r u) dr' 


(4) 


Now using Equation 2 for Gfr/f) and recognizing that the following relationships apply when T' is 
on the scene plane (z' • 0): 


Qr/r') ■ 0 

•^G . 

R « If -Vl 


Equation 4 becomes: 


V(r.») • 2 




scene 


(5a) 

(5b) 

(5c) 

(5d) 

(5e) 

( 6 ) 


wtim E|(?. e) are the electric field componoits on the scene plane which ate tangent to the plane. 
Equation 6 to an expression for the voU^ r^tmse to radiatitni from the scom plane. V( r . e), of 
an antenna with power pattern PqU/F) • (l/A^) A(f/F, v) • A*(f/F. p). The solution as given in 
Equation 6 to In the fiaquency domain. The time donurinrequune to dniredfbr the analysis to be 
performed here and will be derived in the following section. 

5 









III. THE TIME DOMAIN ANTENNA RESPONSE 
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To obtain the time dependent form of the antenna response V(r,t) it is necessary to take the 
Fourier transform of Equation 6. In doing so, it is convenient to introduce the complex analytic 
voltage Vq (r, t) defined as foQows (e.g., Bom and Wolf, 1959): 

ee 

Vc (r, t) ■ 2 J V(f, V ) e^2»r v t<j„ 


The real part of Vg(r, t) is the actual (Le., physically measurable) voltage out of the antenna. 
To make the integration in Equation 7 tractable, the following three approximations will be 


First, it will be assumed that the distance between the image and scene planes is much greater than 
the maxi mum extent of either. The geometry is illustrated in Figure 1. The origin of the coordinate 
syston is in the scene plane at x' ■ y* ■ z' ■ 0 and the position of the antenna is measured from an 
arbitrary point (X^, in the im^ plane by the local coordinates Antenna 

positions are beu^ described with respect to (Xq, Y^, Zq) to accommodate antennas which view the 
scene obliquely with their main beam pointing ak>i« Rq. Zq is the distance between the image plane 
and scene phne and the assumptkm beh^ made here is that (x')2 (y')^ and 

wiaii compared to Z§. This allows the dtetance R between the antenna and an arbitrary point on the 
scene plane to be approximated by: 

ttiilf) • |(Xo-»-»tx**V>(Yo+ny.y')2-*-(Zo + »»2)2p 

• (iix-*’)cos« ♦(i|y-y’)coe(l + i|j«»s7 x-x’)2 + (i|y-y')24-f|2 j 

where cos a -Xo/Ro.cosU - Yq/Rq and cos7 ■ Zq/Ro direction codnes of the vector 

from the or l gto fai the scow ptane to the lefnenM point (Xq, Y 0 , Zg) in the image plane. With this 
approximation and assuming thm IcRq » 1 Equation 6 on be written: 






Second, it will be oanimed that the signal ffom the antenna passes through a filter, h( » ). repre- 
senting the effective bandwidth of the system. This might be an actual fQter placed in the system 
or might r e pr ese nt the effect!^ bandwidth of the detector, integrator or other piece of hardware 
in the system. In either case, it will be assumed that the frequency limitations of the system can be 
described by a filter at the antenna terminals. In addition, it will be assumed that h( v ) is a narrow 
bandpass filter centered abou» frequency p q. That is 


{ H( • Vq) IV. VqK a 

( 10 ) 

0 I V . I > A 

where L/p q<< I. Thb apinoximation means that the receivtav syMem is tuned to receive 
oiersy only in a narrow band about the center frequency p^ With this Mmmption the output 
of the antenna-filter system is V(f, v) H( v • v^). 

Finally, it will be anutned that in the fiequmcy interval passed by the filter, the time and space 
dependence of the electric flehls on the source are separable and that the charactcdstics of the 
antenna dm't clnnge appiecMrly: Thatia.l^?,e) ■ l|(?)B(ir*VQ)and It A(f/r ;e) a 

kg A(F/f ; Vg) whoe ■ 2n pje is the wave number at the center frequency of the filter. 

•• 

Now, with the precadir^ approafanationg the hrtegration in Equatkm ? on done analytically. 
Substitutii^ Equation 9 into Equation 7 and mdcing tin chanfs variaMn I ■ v • one obtafaiB: 


1 


7 
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09 

Vc(?,t) - 4je*i f fj V(f/f i; o) B{ I ) Ha ) exp I [t . r' ) ] j dr'd| 


scene 


where 

V(f/f', j»o) - Et(r') • A(r/r' r) (12) 

and t ■ t • Rg/c is time “retarded” by the propagation delay between the origin in the scene plane 
and the reference point (Xq, Yq, Zq) in the im^ plane. Since H(|) is zero for |{ I > A and 
since A « o the lower limit in Equation 1 1 can be formally extended to infinity. In this case 
the int^ral over t becomes a Fourier transform and one obtains: 


Vctr,t) 


2 w V 


V(f/r V q) B(t •♦/c) • H(t •♦/c) dr' 


where the arterM (*) denotes a convolution and the tilda ( ^ ) over B and H has been removed 
to indicate a Fourier transformed quantity (ie.. the Fourier transform of B(t) is B(t')). 

Equation 13 is the complex analytic form for the response of the anteiun to radiation incident 
from the scene plue. It reduces to standard results in special cues. For example, comdder radiation 
ftom a shigle point source at ?’• ?o . In thia case Vg(r, t) • 2 V(r/f^. ‘ 

exp(>j 2 9p0 Tq) where ■ t « R^c • if (ii,t |^/c is time delayed by the propi^tion time 
from the source at Fq to the antenna ud9(f/ro, no)* 2 l [ko«»‘>^/^ P*X(r/r^ »o)is, 
excqrt for a plm factor, the frequency donnin req^mise of the antenna to a pohet source at ?o 
withpotariationp. When the (Bter Is very narrow [ La.. H(ll« 8(1 )J then V^r,t)« 

2 V(r/fo, no)i(u 0 ) exp(«j2 e which Is a sinusoid at the center ftequency of the Alter 

(Vq) with an amplitude equal to the re^tme of the antrana at fteqoency Uq tiroes the Hwctrum 


I 

{ 

» 

\ 


I 



B 




of the source at v^. In the other extreme, when the filter is very broad (i.e., H(t) * 8(r)) , then 
Vg(f, t) s 2 W(t/Tq, exp(-j2 ir Vq Tq) which again is a sinusoid at the center frequency 

ot the Alter but in this case modulated by the time dependence of the radiation from the source. 


9 
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IV. SPATIAL CORRELATION FUNCTION 

Now consider two antennas at positions r] and r2 in the image plane and suppose that the signals 
Vc(f,t) from each antenna ate multiplied and averaged to form the spatial correlation function 
< V^fi,t)V^r 2 ,t)>. In forming the averages it will be assumed that the spatial and temporal 
variations of the source are independent random precedes and that each random process by itself 
is uncorrelated. That is, it will be assumed that: 

<B(nB*( {')> - f(5 ) 6 ( J - {') (14) 

and that 

<it(r)Et(r')> » ?(f) S (f - ?) (15a) 

S(r ) in the preceding exp’^sion is a dyadic correlation function defined by 




S^(r ) y X 


?(f) 


(15b) 


|^Sxy(f)xy Syy(f)yy J 

where S||(r) is the cross correlation function for the i and j components of the electric field on 
the source. 


Now using Equations 14 and 15 and Equation 1 1 for V^f,t) and doing the integrations over the 
delta frinctions, one obtains: 

OO 

<Vg(fi.t)Vc*(r2,t)> * r JJ dr'J* d I ^A(fi/f vq ) • ?(r') • A*(f2/r,'vo) j 

scene >oo 

( 16 ) 

T({)H2(|)exp|j-^(V» f>] | 

where H^(|)bH(S)H*({)- Equation 16 can be ^mpUfled by letting ^ ^nii 
« r and noting that Ar»-^^Xo*x' ] [Yo-yJ+^^Zo + Slfn i,q '»> (17) 


10 


< 


y 


•f. 


mn 


y 



where 


^ V I * ^ Y ^ 

*1 ^2 V 
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(1 8a) 


‘'x = 


cRo 


^yr^V2 ^ 

c Rq 


^ z\' ^ zn „ 

V. M — i — V 


(18b) 


(18c) 


cRq 


C(*J 1,1)2) 




(18d) 


2 cR(j 

Ar is the difference in time of arrival at the two antennas of radiation from a point on the scene. 
Doing the integration over | (a Fourier transform). Equation 16 becomes 

< Vc(fi,t)Vc*(r2,t)> * ^ J e |*[ A(ri /r , ■ S(r ) • A*(r2/r', *» o)j 


scene 


(19) 


T( At) • H^(At) exp| j2ir ^*'x ^ ‘'y(^o ’ y ) ■*" ’'z^oj | 


where T( A t ) is the temporal correlation function of sources on the scene plane and H- ^ A t ) - 
H( At) * H*(Ar) is the convolution of the temporal response of the filter with itself. 


Notice that if the antennas are identical and close together compared to R^, conditions likely to be 
encountered in remote sensing from space, then Equation 19 can be amplified. In particular. 
i^ssiiming that both antennas are identical with polarization in the $ • direction and so close together 
that their antenna patterns on the scene are essentially the same (Le., MtiI? , ' A(r 2 /r', • 

A(ro/f o) P where r^ is the position vector of the reference point (X^, Z<,) in the im^ 

plane) one has: 




where Pn(r/r') is the normalized power pattern for the antenna and Ag is its effective receiving area 
(Kraus, 1966). Also notice that at microwave frequencies 100 MHz or so is typical of the band- 
width of radiometers but thermal sources such as the earth have a very much broader bandwidth. 
Hence, the assumption, T( Ar) * H^( Ar) ^ H^( Ar) ought to be reasonable for microwave radi- 
ometers viewing the earth. Asuming the antennas are close together and letting T( Ar ) * H^( Ar) = 
H^( Ar) Equation 19 can be written in the final form: 


Vc(r’i.t) Vc*(f2,t)> = 4-^^ e^^°“ f f H^fAr) P„(fo/r") [$ • 5^) • p 1 

scene 

exp|j2ir [(Xo-x')Vx + (Yo-y') «'y+ *'zZo]j 


where $2, = X^/A. is called the anteima solid angle (Kraus, 1966). 
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V. SOLUTION FOR INTENSITY ON THE SCENE 


cr 

Cr 


4 ' ’J 


It is possible with restrictions which are reasonable for microwave sensor systems to invert 
Equation .1 to obtain an expression for p ■ §(r V P> To do so, first introduce the function P (f j “ 
A (r ) where A (T ) « 1 for I? I< b and z«$ro otherwise. PlV) is an idealized antenna pat* 
tem wh'ch neglecu adelobes which are sufficiently far from the main beam. Secondly, assume that 
A V ) a (0) (or r < b. This assumption means that the Alter response is essentially constant 
over the entire Add of view of the antenna. It imposes an upper limit on the bandwidth of the 
measuns' lent. This upper limit depends on the maximum spacing between the anteiuias and varies 
I 'om about 3 GHz * (antenna separation in meters) when to about .3 GHz ^ (antenna 

separation in meters) when Rq » Zq. VITith the Arst of these two approximations the limits of 
integraticn in Equation 21 may formally be extended to infinity and by Imposing the secmid 
approximation the resultant integral takes the form of a Fourier transform on v and t' y. 

With these two approximations Equation 21 is a Fourier transform which can be inverted to 
obtain an expression for p*¥(r) ’ One obtains 


P ^ ^ ° Jf * e *^220 < Vc(f|.t)V;*(;2.0> 

-OO 

exp j rj2rr [( "x * "y ] ) ***'x 
where Cq * 008^7 K^CO)] . 

Eadi of the terms L :he integrand in Equation 22 (le. ('x# I'y* 

from the po*>Jons ri and ?2 of the antenims. Consequently, if < V^q.t) V^*(f2,t) > was known 
as a ftw .tion of position, the integrand in Equation 22 would be determined and the integration 
coi..«« be peiformed (e^. mmerically). In (bet, < V^ri,t) VQ*(r2,t) > can be measured udng a 
ooherait detecuw sudt n diown in Figore 2. Referring to P^ure 2, V|^2(^) cosfeut ♦ 9j 2^ • 
ReVc(r| ^ t) ate the real (Le. phyaiGdlymeanitd>te)volti^attheantainatennimitoandthe 
dedtfi conation Ametion is <V^fi,t)Vc*(r2 ,t)>« <V|(t)V2(t) [ cos(9| > 02^ 

13 



jsin(0| * 02 > 3^* obtain the correlation function the real voltages from each antenna are 
initially coherently mixed with the local oscillator signal to shift them to a carrier c«;|p. Then 
the signals are multiplied together and low pass filtered to obtain V j(t) V 2 (t) cos(0| - 02> which is 
the real part of V(.(f- , t) V^*(r 2 , t) and also multiplied after one signal has been shifted in phase 
by ir/2 radians and then filtered to obtain the imr^ary part of the correlation function. Assum- 
ing that Vf^r,t) is a stationary, etgodic random process, the ensemble average (denoted by the 
pointed brackets <> ) can be obtained as a time average. Thus, given a measurement in the image 
plane, the v and $2 are determined by the position of the antennas and < V^r^, t) Vc*(r2, t) > is 
determined by avert^g the signal from a coherent detector such as shown in Figure 2. Given these 
quantities the integral on the right side in Equation 22 can be evaluated numerically to determine 
the intensity of radiation on the scene plane. 

Figure 3 illustrates how this procedure m^t be implemented to do remote sensing from space. 
Imagine two identical anteniuis, one spiraling about the other, and suppose that the signal from 
each is received at a common point and processed as indicated in Figure 2. At discrete points indi- 
cated by the dashes «n Figure 3, the average < V^fi,t) V^*(r 2 ,t)> is recorded. Since each point 
corresponds to a different baseline, aiu therefore to different frequencies ^ y» ^ S 3 rtem is 
in effect mapping the correlation function < V(^ri, t) Vc*(^2» t)> as a function of v ^ and v y. 
If the coordinates ^ x> ^ y> known at eadi point, then all the terms in Equation 22 are 

determined. After all the measurements have been made the integral can be evaluated numerically 
to determine the intensity of radiation on the scene plane. The accuracy with which this can be 
accomplished will depend on the range of Vy over which measurements have been made 
(the diameter of the spiral in Figure 3) and the spacing between measurements, but it does not 
depend on the size of the antenna used for this measumnent. The antenna must be diosen large 
enough so that the paraxial ray approxhnation (Equation 8) is valid, but otherwise the antenna 
size does not affect the resolution of the measuntment A large antenna reduces the field of view 
(area imaged) but does not determine the resohithm. 


14 


The scheme presented in Figure 3 is conceptually what must be done to obtain the intensity on 
the scene plane using aperture synthesis. The essential requirement is to make measurements at 
different inteferometer baselines and the important flexibility is that these measurements do not 
have to be made at the same time. Thus, one can move a single antenna pair sequentially from one 
baseline to the next as in Figure 3, or one could have placed antennas at aU the dashes and obtained 
measurements at all baselines at once. There are many possibilities between these extremes. For 
example, antennas uniformly ^aced along two orthogonal axes, perhaps in the form of a +, an 
arrangement called a Mills cross, can measure all the baselines in a rectangular grid of the same 
dimensions as the cross. Additional baselines would be possible if the Mills cross was slowly 
rotating. It is also possible to keep the anteima positions fbced and change the baseline by chang- 
ing the frequency of the measurement (Equation 18). Aperture synthesis using this approach has 
been su^ested for remote sensing from space by Schanda (1976), but it requires measurements 
over a large frequency range to obtain good resolution. Other variations of the basic concept have 
ateo been proposed (Mefnik, 1972; Wiley, private communication). 


15 
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VI. AN EXAMPLE 

To Ulustrate the procedure consider i system in which measurements are uniformly spaced on a 
Cartesian grid. That is, assume that the antennas in the image plane can be at positions (v = n d ; 
fly * «ydy ; Vz • 0) where n^ and ny are integers. Then assuming that | n^ y I < y and letting 
^ 0 *^ 0 *^ fo*" convenience one can write Equation 22 in the following form: 


N, 


N. 


Ipp(x'.y ) - J J V(nxAx. nyAy)ej2'[nx^xX+ OyAyyl (23) 


"x*-Nx ny»-Ny 


»*««V(r,.,y).c„e-»oO<vrff,.,)V,-<Fj,,)> Comi»ra,Equ» 

tioiis 22 and 23, one sees that tpp(xjr ) is the intensity on the scene plane with polamations ii 
Obtained with this particular sampling scheme and V(v^, Vy) is the Fourier traiwform of the 
actual intensity on the surface multipUed by the normalized antenna pattern. That is, 

F -I [ V( V g, vy)] « p(f') . sCr ) . - P(?') ipp(x', y') where Ipp(x’,y ') is the true intensity 

on the scene plane with polarization in the ^ direction and P» denotes an inverse Fourier trans- 
form. The significance of Equation 23 can be made more apparent by introducing delta functions 
* ”x,y ^,y)' Doing so one obtains: 

Wv)- ^2//o»,J-yV(p„p,,c-i2'[«,xH.-yy]„ 

"x,y 


» Fkr^d «»xdi» yV(t>^.Py)[^ Ay J 8(v X* nxAx) 8 ( «'y nyAy)J e ^2 * [•'x* +»yy] 

"x.y 

- r Vr ) F-l [ V(,,, ,j,)J . D(V,y ) 


The (•) d«iow . coimihitton and D(x^ y') i, tte Fomiw timsfom of the sum 

Ax *yi»<'x-"xAx)»(»y-iiy4j,)whWiceiibewrttteiitathefollowhiifofin: 

^.y 
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D(x'.y') 
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» h(x)g(y’) 


sin[ir x'Ax(2Nx+ D] 
* sin[irx'Ax] 



sin [ IT y'Ay(2Ny + 1)] 

sinC^y'^y] 


( 25 ) 


Notice that Equation 24 has the fonn associated with the response of an antenna with effective 
(power) pattern D(x',y') (e.$., Kraus, 1966; Collin and Zucker, 1969). In fact, h(x') and g(y') are 
just the array factors encountered in antenna engineering for a uniformly spaced, linear array 
(Jordan and Balmain, 1968). Thus, in this example the measurement described in Equation 22 
yields an image which has the same resolution obtained with a planar array of (2Nx,y * 1) 
elements on a side spaced dx,y meters apart The important point to notice, however, is that this 
resolution is determined by the pattern of measurements taken in the image plane, and not by 
the antenna used in the measurements (Le. not by Pn(r/r') ). This means that measurements with 
high resolution can be made even with small anteimas by making measurements over a sufficiently 
large area in the image plane. 

Finally, notice that in the limit \ y N^ y -*■ «o, the functions h(x') and g(y') become 
delta functions at x' * y' « 0. In this casetpp(x', y ) is equal to the actual intensity on the surface, 
Ipp( x', y' ). In this case the measurement is perfect. In the more general case h( x' ) and g( y' ) 
have peaks of finite width (Nx,y^x,y^*^ y ) h a blurred image of the true intensity, 

Ipp(x', y'). The width of the functions h(x’) and g(y‘) is a direct measure of the blurring and 
can be used as a figure of merit for the resolution of the system. Thus, using the width to the 
first zero to define the resolution of the synthesized antenna aperjure, one has: 

^x,y • [Nx,yAx,y]*^ 

• ^ ^ 

Nx.y<*x,y •'o 

Although the antennas used in the measurement don’t determine the resolution, they do determine 
the field of view (FOV) of the synthesized antenna aperture. Assuming measurements made with 
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antennas with din.ensions 0,^ y in the x- and y-directions, respectively, the portion of the scene 
plane from which energy is received is approximately; 


FOVv 


n ^o 
*^x,y 


c 

«VDr,y 


'x.y ■ Dv„ *'0 i?TU7v (27) 

where is the distance from the antenna to the scene plane. However, several restrictions have 
been imposed on this Held of view in the course of developing the theory. Pint, the field of view 
must be much smaller than the distance between the image plane and scene plane to satisfy the 
sagittal approximation (Equation 8) made in the derivation of the antenna response (Equation 9). 
Using Equation 27 this restriction may be stated as: 

c 


^x,y 0 


« 1 


(28) 


(The sagittal approximation also limits the size of the array possible in the image plane, but con* 
sideling that Rq is large for sensors in space, tha is not likely to be a significant restriction.) Second, 
the array functions h(x' ) and g(y') have grating lobes which affect the field of view. These grating 
lobes occur every lMx,y nieteis; hence, for unambiguous measurements the field of view must 
be restricted to be less than this distance: 


FOVx^y < "2^ • Nx,y fix,y 


(29) 


Comparing Equations 27 and 29, one finds that dx,y <Dx,y to avoid grating lobes. 


Combining these results, one obtains the following parametem for the ima^ formed with this 
measurement scheme. Assuming measurements from an altitude. R^, at a frequency, with 
antennas of size, D^^y on a side, and assuming 2Nj^ y 1 measurements uniformly spaced d^^y 

meters apart on each side, one obtains: 

Resohititm 


P« 




(Nx,y dx,y) *'o 

Field of View 


FOV 


x.y 


cRc 


(30a) 


(30b) 


18 



For no Grating Lobes 



VII. CONCLUSIONS 






Notice that in the preceding example. Equations 30a-c are just the parameters associated with the 
antenna beam of a linear array with uniform element spacing d^ y and N^ y elements on a side. 

That is, the image synthesized with the uniform sampling scheme described in the example will be 
the same as one formed b]^scanning with a linear array which occupies the same region in the image 
plane as was used to obtain the synthesized image (Equation 22). However, the unique advantage 
of the synthesis procedure is that all of the antennas do not need to be present for each measure- 
ment. In principle, the image could be formed with just two antennas which are moved to appro- 
priate positions in the image plane. 

Also notice the special role played by system bandwidth in the synthesis scheme. The field of view 
obtained with aperture synthesis and with a linear array with an equivalent resolution are determined 
by the elemental antennas used in the measurements; however, in the case of aperture synthesis the 
field of view is also affected by the bandwidth of the measuring system. For example, consider an 
ideal filter with response H^(^r) • H^(0) for £ir < and 0 otherwise. Radiation from a point sym* 
metrically located below a pair of antennas (Le. equally distant from each) arrivn at each antenna 
at the same time (le. * 0); however, as one moves away from this position the time delay 
increases. When the delay exceeds Tq radiation from this position no longer is passed by the filter. 
Cmsequently, the filter has the effect of restricting the field of view. To a first approximation 
(Xq • Yq • 0 and much larger than the field of view) this distance is cR^/d^ y A where 
A is the system bandwidth (A a r and d^y is the distance between antennas. In the 
analysis presented above it was assumed that H^( Ar) ■ H^(0) for points within the field of view 
determined by the antennas. That is, the field of view set by tlm filter was chosen to be much 
greater than that set by the antomas thermdves. However, it hi also posribto to use the filter itself 
to control the field of view. The difficulty fe that At, and therefne the field of view, depend (m 
the spacbig between the antennas. Consequently, to keep the field of view conttant. the filter 
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must be changed for each measurement. This may not be so impractical as it seems if, for example, 
the data from each antenna was resm’ded digitaUy with a fixed bandwidth and then processed 
later to obtain the correlation fiuiction. The filtering in this case could be done numerically as part 
of the processing. Such a procedure would have very great flexibility in choosing the field of view, 
but it requires very large band width during the initial recording process. 


Another point to consider is the sensitivity of the synthesized aperture (Le., the accuracy with 
which the source can be determined with a given averaging time per measurement). To give an idea 
of the sensitivity achieved with aperture synthesis a comparison with the sensitivity of a linear array 
has been nmde. This is presented in the Appendix. It is shown in the Appendix that for point 
sources the real aperture is more sensitive because its narrow beam restricts the background against 
which the source must be contrasted, whereu in the case of aperture synthesis the antermas view 
the entire field of view on each measurement. However, for extended sources such as the earth, 
which fill the field of view, the real aperture and synthetic apertures receive the same power density 
on each measurement and therefore are equally sensitive. Hence, for remote mnsing of the earth 
from space, aperture synthe^ can provide sensitivity comparable to a real aperture. 


RnaDy, coorider the effects of antenna portion eiron on the imc«e obtained with the synthe;iis 

sdieme. An error Ai|| in antenna separation along the i>th coordinate is equivalent to an error 

AV| in spatial fieqoency (Equation 18). In ordm for this error to be have a smaD effect on the 

/AdiX Xi 

image one requires (Equation 22) that X| •( ^ where is the wavelength 

at vidiidi the memuament to made and X^ 9 |(Xo ^« ). ^ V*)* * However, to satWy 

the sagittal ^pro a d ma t hm (Equation 7) it hn bean aaniieed that 1 mid y'/Ro « 1. 

Hmce, the c ondi t to ntiMtpt^ttoienontewesmaH effect on the hnagawn be written 
« 1 whoa X{ e |XoiYo»Zo}*NociathttXo/Ro*^^ 


with mdto. He nce, the synthetk aperture on be made hontotive to errors in the position parallel 
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to the image plane ^Vy)^V restricting the antennas to point near nadir. On the 

other hand, except in extreme cases Zq ^ Rq. Consequently, em>n in 'he position perpendicular 
to the image plane are particularly critical and need to be much less than a wavelength. Similar 
comments pertain to the linear array. However, notice that the requirement in the synthesis pro- 
cedure is onl)' for knowledge of position, and that the position itself is not particularly important. 
If the position is known, even if it is many wavelengths from the desired position, then appropri- 
ate correctk>ns can be made in processing .he data by adjusting the phase terms in Equation 2?.. 



VIII. SUMMARY 


A theoretical basis has been presented for synthesizing antenna apertures for microwave radiometers 
in space. The technique consists of measuring the product of the signal from pairs of antennas 
at many different antenna spacings, thereby mapping the correlation function of the antenna 
voltage. The intensity of radiation at the source can be obtained from this correlatior function by 
means of a Fourier transformation. A variety of procedures are possible to measure the correla* 
tion function, including a pair of antennas with variable baseline, an array of anteimas which simul* 
taneously make m . . surements at many baselines, measurements with fixed spacing but at several 
frequencies, and combinations of these. For mapping the earth from space all of these schemes 
must limit the field of view in order to guarantee a Fourier transform relationship between the 
intensity on the scene and the correlation frmction. The resolution obtained with these schemes 
is amilar to that of a real linear array of the same dimensions in the image plane. The mqjor 
differences being that the synthesized image is obtained by sigrial processing after all measurements 
are made and that the field of view of the synthesized aperture must be restricted to meet the 
paraxial ray approximatiotL Ihe field of view can be limited by the actual anteruias employed 
in the measurement or by controlling the effective bandwidth of the measurement. The latter 
requires large bandwidth but is potentially the most flexible system. 
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APPENDIX 


Although the synthetic apertute has the same resolution as a linear array of comparable size, the 
images are formed in a much different manner with each system. The linear array maps by forming 
individual beams which i-eceive energy from each resolution cell in the map (i.e., one beam per 
pixel). In contrast, the synthetic aperture makes all its measurements with a small antenna which 
receives energy from the entire field of view during each measurement. This raises the question of 
the relative sensitivity of the two approaches: that is, given an available integration time per 
measurement, are there differences in the temperature resolution of the images formed with 
aperture synthesis as compared to that obtained with a real array. This question is addressed in the 
paragraphs to follow where it b shown that for point sources the real antenna is more sensitive but 
that for extended sources which fill the field of view, aperture synthesis and the real array are 
equally sensitive. 

• 

For thermal sources the power available from an antenna with normalized power pattern, Pq (9, 0) 
can be expressed in terms of the equivalent black body temperature of the sources in the foUowii^ 
form (e.g., Kraus, 1966; this can also be derived from Equation 9 in the text): 

T(9,0)Pn(9,0)sin 9d9d0 (Al) 

where 3 is caOed the antenna beam solid an^e and equals (directivity/4 ir )*^, is a constant 
called the anterma efficieiicy to account for losses in the antenna system, and T 3 is the noise temp* 
erature of a resistor which has the same available power as the antenna. Assuming for simplicity 
that the temperature is constant over the source, T(9, 0) et T^, the above Amplifies to: 

<«) 


2S 



where 
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“s* JJ 


,0) sin 0 d9 d0 

'source 

The integration time r required to measure with an RMS error (i.e., variance) is deter- 
mined by the system temperature and bandwidth A in the following form: 


(A3) 


e T 


6Ta 


sys 


where e is a constant which depends on the type of receiver ( e s: .7 for a correlating receiver 
and £ SI 2 for a Dicke receiver). Combining Equations A2 and A4, one obtains the following 
expression for t : 

e X 


(A4) 


/ e Tjyj \ / Q a \ ^ 

\ A 


(A5) 


n&To 

It is desired to compare the time required to form an image using a large real antenna and a syn- 
thetic aperture with the same resolution. To form the synthetic aperture assume that measurements 
are made uniformly on a cartesian grid with y 1 measurements on a side and in order to 
obtain the maximum unambiguous field of view assume that d^ y » y where y are the 
dimensions of the anterma used in the synthesis. Let S2g be the solid angle of this antenna and let 
be the solid angle of the equivalent real aperture. The real aperture is y d^ y meters on a 
^de and neglecting possible differences in side lobe levels definitions, 

the time required per measurement with the real antenna is: 


/^A^sysV/nA Y ^ 


,»a‘\ 

and the time required for each measurement in forming the synthetic .«.«iture is: 

a 


(A6) 


^a ® 


/ *a^sys \ ^ ^ \ -i_ 


(A7) 
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where and are the beam solid angles subtended by the source in the case of the large 
antenna and the antenna used in aperture synthesis, respectively. 


The total time required for each system to obtain an image is obtained by multiplying by the 
required number of measurements. Assuming that these are the same in each case (actually, fewer 
measurements are needed in synthesis schemes which remove redundant measurements), and 
assuming that each system has the same bandwidth and approximately the same system temperature 
and efficiency, ore has: 

^sa 
^SA 

where it has been assumed that 2N^ • 2Ny « N for simplicity, and the relationship £2^ » £2 
has been used. 

The ratio in Equation A8 may have several different values depending on the size of the source. For 
example, when the source is smaUer than the main beam of either antenna ^ SA ‘ ^ sa> 
source fills the beam of the real anteima but is smaller than the beam of the antenna used in 
synthesis £2§y^ = £2^^; and when the source is larger than the beam of either anteniu £2§^ » 

£2 A and £2 sii a £2 a. Thus, one obtains the following possibilities: 








source smaUer than either 
antenna beam 

source fills beam of real 
antenna but less than FOV 

source fills beam of both antennas 


Clearly, in the case of anall (e.^, point) sources mapping by scanning with a single large antenna is a 

. N 

more effidbnt use of available energy than aperture syntheds and requires 1^ integration tinm per 
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measurement to obtain desired accuracy. However, when the source completely fills the beam of 
the antenna used ) i synthesizing an aperture, then aperture synthesis and scanning with a large 
antenna yield comparable sensitivity. This is the situation likely to be encountered in remote sensing 
of the earth’s surface. 
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calculations. 
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V^(t) cos M + 


V2<t) COS (cot + ^2) 


ER 


LOCAL 

OSCILLATOR 


LOW PASS 
RLTER 


LOW PASS 
RLTER 


Vi(t) COS (<t>|Ft + «»^) 


V2(t) cos (Wipt + <P2) 


LOW PASS 
FILTER 


PHASE SHIFT 
90 * 




LOW PASS 



RLTER 

1 

V, V 2 «ln( 9 , -9*) 

T 

V| V2 COS ( 9 ^ — ^2) 

F^ttf»2. R«oriver for m^urii^ tbe cofielation Auiction. 


30 











